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Abstract 

The Schrodinger equation is investigated in the Euchdean Taub- 
NUT geometry. The bound states are degenerate and an extra de- 
generacy is due to the conserved Runge-Lenz vector. The existence 
of the extra conserved quantities, quadratic in four-velocities imphes 
the possibihty of separating variables in two different coordinate sys- 
tems. The eigenvalues and the eigenvectors are given in both cases in 
explicit, closed form. 
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1 Introduction 



The discovery of the self-dual instanton solutions to Euclidean Yang-Mills 
theory suggested the possibility that analogous solutions to the Euchdean 
Einstein equations might be important in quantum gravity. One example 
of a metric which satisfies the Euclidean Einstein equations with self-dual 
Riemann tensor is the self-dual Taub-NUT metric [Q. In this case Einstein's 
equations are satisfied with zero cosmological constant and the manifold is 
M.^ with a boundary which is a twisted three-sphere possessing a distorted 
metric. 

On the other hand much attention has been paid to the Euclidean Taub- 
NUT metric because it is involved in many modern studies in physics. The 
Kaluza-Klein monopole of Gross and Perry |^ and of Sorkin [§] was obtained 
by embedding the Taub-NUT gravitational instanton into five-dimensional 
Kaluza-Klein theory. Remarkably, the same object has re-emerged in the 
study of monopole scattering. In the long-distance limit, neglecting radiation, 
the relative motion of two monopoles is described by the geodesies of this 
space |, §. 

It is well known that space-time isometries give rise to constants of mo- 
tion along geodesies. However it is worth to mention that not all conserved 
quantities arose from isometries of the manifolds and associated Killing vec- 
tor fields. Other integrals of motion are related to "hidden" symmetries of 
the manifold, which manifest themselves as tensor of rank n > 1, satisfying 
a generalized Killing equation, namely V(^-ft',yi...i,„) = p. They are usually 
referred to as Stackel-Killing tensors. 

An illustration of the existence of extra conserved quantities is provided 
by the Taub-NUT geometry. Indeed, for the geodesic motions in the Taub- 
NUT space, there is a conserved vector, analogous to the Runge-Lenz vector 
of the Kepler type problem, whose existence is rather surprising in view of 
the complexity of the equations of motion |]^, ^, ^ |10[. The Runge-Lenz 
vector is quadratic in four-velocities and its components are Stackel-Killing 
tensors. 

The Taub-NUT space is also of mathematical interest. In the Taub-NUT 
geometry there are known to exist four Killing- Yano tensors. The Killing- 
Yano tensor [11] here is a 2- form, f^^ = f[^u], which satisfies the Penrose- 
Floyd equation V {^fu)x = 0. Three of these Killing- Yano tensors are rather 
special: they are covariantly constant, mutually anti-commuting and square 
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the minus unity. Thus they are complex structure reahzing the quaternionic 
algebra and the Taub-NUT manifold is hyper-Kahler [^. The fourth Killing- 
Yano tensor has a non- vanishing field strength, is not trivial and leads to new 
constants of motion. More precisely, the symmetrized product of the fourth 
Killing- Yano tensor with the previous three is connected with the Stackel- 



Killing tensors, i.e. the components of the Runge-Lenz vector p, [T^, |T3|, |14 

The aim of this paper is to study the Schrodinger equation in Taub- 
NUT space. The conserved Runge-Lenz vector implies an extra degeneracy 
of the bound states. It is remarkably the fact that there appears to be a 
close relation between the existence of extra conserved quantities and the 
possibility of separating variables in two different coordinate systems. In 
this way one obtains two different orthonormal energy bases, called here the 
central and axial bases, for which the degeneracy of the energy levels is well 
defined by the eigenvalues of suitable sets of commuting operators. 

The plan of the paper is as follows. In Section 2 we summarize the 
relevant features of the Euclidean Taub-NUT geometry. In Section 3 we 
analyze the Schrodinger equation in this space looking for the complete sets 
of commuting operators that define the mentioned bases. In the next two 
sections the Schrodinger equation is investigated in spherical and parabolic 
coordinates. We show that the central modes can be solved in spherical 
coordinates while for the axial modes we need to use parabolic coordinates. 
In both cases we obtain sets of regular modes as well as specific sets of 
irregular modes arising from the special boundary conditions of the Taub- 
NUT geometry. Our concluding remarks are presented in the last section. 
Finally, the Appendix is devoted to S0{3) ® U{1) harmonics. 



2 Taub-NUT geometry 



The Euclidean Taub-NUT manifold is a particular case of static Euclidean 
Kaluza-Klein space-time whose metric, in a static chart of coordinates, t and 
x^, (/i, z/, ... = 1, 2, 3, 5), is given by the line element 

ds"^ = -dt^ + g^^dx^dx" = -dt^ + ^df + V{dx^ + Aidx'f , (1) 

where 
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is the usual Euclidean 3-dimensional line element involving the Cartesian 
space coordinates ... = 1,2,3) which cover the domain D. The other 

coordinates are the time, t, and the Cartesian Kaluza-Klein extra-coordinate, 
G -D5. We suppose that V and Ai are static functions depending only on 
X. It is clear that in the absence of the scalar potential we have 

div A = , TotA = B. (3) 
/^From (|l|) we obtain the covariant components of the metric tensor 

gij = ^Sij + VAiAj , = V Ai , g^^ = V , (4) 

and ^ 

g = det{g^^) = — . (5) 



The corresponding contravariant components are 

g'^ = V5,,, g'' = -VA,, g^5 = ± + vT. (6) 



1 -.-2 
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The Euclidean Taub-NUT manifold has the line element of the form (|^) 
with 

^-^ = 1 + ^, ^1 = --^, ^2 = ^^, ^3 = 0, (7) 

where r = | x | and /i is a real number. Hereby, it results the magnetic field 
with central symmetry 

fi = /ij. (8) 

In fact, the geometry defined by (|T]) and (|^ has the global symmetry of the 
group S0{3) (S> f^5(l) ® Tt{l). This means that the line element is invariant 
under global rotations of the Cartesian space coordinates and x^ and t trans- 
lations of the Abelian groups U^i^l) and Tt{l) respectively. The generators 
of these groups are just the five Killing vectors corresponding to this global 
symmetry. Moreover, there is a three-component Killing tensor kfiu which 

satisfies V(o- k^v)= (where V is the covariant derivative and ( ) denotes 
the symmetrisation). 
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3 Observables 



The quantum mechanics in the Taub-NUT geometry ||T0[ is based on the 
Schrodinger equation 

Hip = idtip (9) 
involving the Hamiltonian operator 

H = -^V,g''''V,, (10) 

written in natural units with h = c = 1. In any static chart, Eq.(^ has the 
particular solutions 

Mx)=UEix,x')e-'^\ (11) 
where Ue are energy eigenf unctions, 

HUe = EUe . (12) 

The solutions (|Tl|) may be either square integrable functions or tempered 
distributions on D x D5. In both cases they must be orthonormal (in usual 
or generalized sense) with respect to the scalar product 

{ip,ip')=f d'x^r^'. (13) 

We denote by £^ the Hilbert space of the square integrable functions with 
respect to this scalar product. 

The main operators on can be introduced by using the geometric 
quantization |1^. In this way, one obtains the non-hermitian momentum 
operators in the coordinate representation, 

= -t{d, - Ad,) , P, = -td,, (14) 
which obey the commutation rules 

[P„ P,] = te,jkBkP5 , [Pi, n] = . (15) 
With their help we can write 
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Hereby we see that the U^ll) generator, P5, is conserved since it commutes 
with H. This is natural since P5 is up to the factor —i just a KiUing vector. 
In the following it is convenient to replace it by Q = —f^P5- Furthermore, 
one can verify that other three Killing vectors are the components of the 
angular momentum operator 

2 = xxP+-Q. (17) 

r 

These are conserved and satisfy the canonical commutation rules. Moreover, 
their commutators with the coordinates and the momentum operators are 
the usual ones. Other three important conserved operators are those defined 



by the three-component Killing tensor k^u as |T^, |T0 



^ 1 -^fj-i^ 1 ^ ^ ^\ X ( 1 o\ , . 

K=--Vf,k \/, = -\^Px L- Lx pj-n-iH-—Q^\ . (18) 

Thus we obtain a vector operator which play the same role as the Runge-Lenz 
vector in the usual quantum mechanical Kepler problem. This transforms as 
a vector under space rotations since its components obey 

[Li, Kj] = ieijkKk . (19) 

Now we have to chose the suitable complete sets of commuting operators 

which should define usual or generalized bases for the Hilbert space C^. Here 

two options are interesting, namely (I) that of the central basis given by 
2 

the set {H, L , L^, Q} or (II) the axial basis formed by the common eigen- 
functions of the set {H, K3, L3, Q}. In both cases we must work in suitable 
coordinate systems which should allow us to separate the variables in all the 
differential equations corresponding to the eigenvalue problems. 
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4 Central modes 



In order to study the central modes it is convenient to choose the local 
chart with spherical coordinates, r, 9, 0, commonly related to the Cartesian 
ones and the new coordinate x defined as 



1 

X = — arctan— . (20) 

Here the Taub-NUT line element is 

ds^ = -de + ^{dr^ + r^de^ + sin^ 6 d(f) + /xV(rfx + cos 6 d(j)f , (21) 
since 

A^ = Ag = 0, A^ = cos 6) . (22) 

Moreover, we consider that in this chart r & D,. = {r\V{r) > 0} (i.e., r > 
if/i>Oorr> if/i<0), the angular coordinates 6, cover the sphere 

while X^D^ = [0,47r). 

The main observables are 



Q = -id^ (23) 



and the components of the angular momentum ([T7|) in the canonical basis 
(with L± = Li ± iL2), 

L3 = -id^. (24) 



L± = e 
The 5*0(3) Casimir operator 



±dg + i{cot6d^ Att dy 

sm 6 



(25) 



= L+L_ + Ll- Ls (26) 
^ de{smede)--^{dl + dl-2cosed^d^) 



sine""' sin'e'"^ ' ^ 
is related to the momentum operators through 

2. 1-2 1 

-Or + - 



P =dr^ + -dr-^L (27) 



Looking for the common eigenfunctions of the set (I), we observe that the 
angular momentum as well as Q act only upon functions of 9, and x which 
must be square integrable since the domain S"^ x is compact. These form 
the Hilbert space £^(5*^ x where the complete set of commuting operators 

2 

{L , L^, Q} determines a suitable basis. Their common eigenfunctions, F^^^, 
which satisfy the eigenvalue problems 

2 

^ ^ l,m 

q 



l{l + l)Y,l 



(28) 
(29) 
(30) 



and the orthonormalization condition 



yl yi 

-' Lmi -' I'm' 



4lT 





Sl,l'Sm,m'Sq^q' , 



(31) 



will be called 5*0(3) ^U{1) harmonics. Notice that the boundary conditions 
on X require / and m to be integer numbers while q = 0, ±1/2, ±1, ... 
[|lO|l . The form of these harmonics is given in Appendix. 

Thus the problem of the angular motion is completely solved and it is 
clear that the angular coordinates can be separated from the radial one if we 
take the common eigenfunctions of the set (I) of the form 



E,Lm 



'r,0,(f),x) 



(32) 



Then from Eq.(|T^) we obtain the familiar radial equation. 



where we have denoted 
a = 



d^ l{l + 1) 



2/i 



Rlmir) 



PRlmir) 



(33) 



(34) 



The radial wave functions must be orthonormal with respect to the radial 
scalar product 



^E,h ^E'l 



f dr(l + ^)R%,{rrR^^,,{r). (35) 
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This is a Kepler-like problem similar to that of the non-relativistic quan- 
tum mechanics. The general solution of Eg . (|33|) can be written in terms of 
the confluent hypergeometric function as 

R%^{r)=Np^e-PF{a,2s,p), (36) 

where 

I Q; 

P = 2r-\/-/3 , a = s- > (37) 

and is the normalization constant. The parameter s is a solution of the 
equation s{s — 1) = /(/ + 1). The modes with s = / + 1 are similar to those of 
the usual non-relativistic case. For this reason we say that these are regular 
modes. We shall see that in the Taub-NUT geometry there is a conjecture 
in which irregular modes with s = —I are also allowed. In the following we 
shall briefly discuss the corresponding energy spectra assuming that E > 0. 

Let us consider first the regular modes (s = / + 1). It is easy to show that 
for /i > there is only a continuous energy spectrum, covering the domain 
E > the levels of which are infinite degenerated as it result from the 

selection rules (§^). Notice that these rules are in accordance with the behav- 
ior of the classical trajectories which are open and cannot reach the center. 
In the case of /i < we have the same continuous energy spectrum with 
infinite degenerate levels and, in addition, a discrete spectrum in the domain 
[0, where a > and /3 < 0. Here we can impose the quantization 

condition a = —Ur, = 0, 1, 2, ... which gives 

= 2n (38) 



where n = n,, + / + 1 is the main quantum number of the regular modes. 
Hereby we find the energy levels 



1 

En = ^ 



nJ n? — q"^ — {n — q 



(39) 



for all n > |g| > [|T^. These levels are finite degenerated since for given q 
and n the angular quantum numbers can take all the values selected by the 
conditions — 1<Z<?7, — 1 and (p5|). Thus one obtains a countable energy 
spectrum with the property 

lim En = ^. (40) 
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The irregular modes arise only as discrete modes when /z < since the 
radial functions (|36|) remain square integrable on = oo) for s = — / if 
a = —rir and |g| + / < n,, < 2/. The main quantum number of these modes, 
n = Ur — I, obeys |g| < n < / and gives the same quantization rule (^) 
and energy levels (|39|) . Therefore, the discrete spectrum of irregular modes 
coincides to that of the regular ones, spanning just the domain of energies 
where the classical motion has closed trajectories. The difference is that in 
the case of irregular modes the energy levels are infinite degenerated (since 
/ > n) and there is no continuous energy spectrum. 



5 Axial modes 



Let us consider now the axial modes given by the eigenvalue problem of 
the set (II). This can be solved in a local chart with parabolic coordinates, 
^ G D^, 1] G and G [0, 2tt), defined by 



X 



^ - '^r]cos(j) , x^ = ^J^r]sm(j) , = ^{^ - r]) , (41) 
such that 

r = i(e + r/). (42) 

Here the domains and corresponding to are either ^, 77 > for 
/i > or ^, ?7 > |/i| if /i < 0. The line element reads 

ds^ = -dt^ + ^df + /iV (dx + 1^ dA (43) 



1 



where 



dl' = ^^de + ^^dv' + er/#^ . (44) 

Following the same procedure like in the case of the central modes we 
calculate how look the operators of the set (II) in these coordinates. It is not 
difficult to show that 

H=^Q'' (X + Y) (45) 
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and 



[(r^ + /i)X - (e + fx)Y] 



where we have denoted 
X = 

Y = 



The whole eigenvalue problem of axial modes is 



E,K,m 5 
? 

E,K,m ' 



E,K,m 1 



These equations can be solved in parabolic coordinates if we put 



(46) 

(47) 
(48) 

(49) 
(50) 
(51) 
(52) 

(53) 



where is the normalization constant calculated with the help of the scalar 
product (O) rewritten in parabolic coordinates as 



(54) 



{U,U') = - f dif dr^i^ + + 2f,) 
4 JDf Jd„ 



X 



An 



dx U*i^, V, 0, x)U'i^, r/, (j), x) ■ 



By taking into account that the operators of Eqs. (^) and {pU\j have the form 
( ^5]) and ( |i6| ) respectively, we find that / and h satisfy 

(m + g)^ 



dm) - 



4e 
[m — qY 
Arj 



a 3 ^ 

+ T + T^ 
4 4 

a 6 
+ T + 



^E,K,miv) 



2 /i?,K,?Tl(0 



(55) 



"2 hE,K,mi^) , (56) 



where a and P are the functions of E and q defined by Eqs. (|3^) . The last 
step is to solve these equations. Defining the new variables 



X 



(57) 
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we obtain the general form of the solutions, 



fk.,m{0 = a;^^e-^/2F(ai,2si + l, x) , (58) 

h''E,.,miv) = y''e-y/'F{a2.2s, + l,y) , (59) 

where 

1 2k — a ,„ . 

a2 = ^2 + --^, (61) 

depend on the parameters Si = ±(|m| + q)/2 and S2 = ±(|m| — q)/2. 

We define the regular axial modes for si = (|m| + q)/2 > —1/2 and 
S2 = {\Tn\ — q)/2> —1/2 which means that |m| > |g| — 1. These correspond 
to the regular central modes we have discussed in the previous section. The 
discrete energy spectrum appears in the same conditions (a > and /3 < 0) 
but now the quantization rules are 

ai = -rii , a2 = -n2-, ni, n2 = 0, 1, 2, ... . (62) 

Hereby we recover the quantization rule (^) with the main quantum number 
of the regular modes n = ni+n2 + \m\ + l and the same formula of the energy 



levels (pQf ). Moreover, we find that 

^ = (^"-^ -^i-q) ■ (63) 

Thus the regular discrete axial modes are labeled by the quantum numbers 
ni, 71-2, Tn and q. Notice that a and (3 depend only on q and n as it results 



from Eqs.dl) and (H) 



The irregular discrete modes can be derived for yU < in the same manner 
like in the case of the central modes, by looking for the conditions in which 
the wave functions (|^) remain square integrable on the domains C,,ti > 
for Si,S2 < 0. We observe that this happens only when ai = —rii, a2 = 
—n2, Si = —{\m\ + q)/2 and S2 = — q)/2 such that rii < \m\ + g — 1 

and n2 < |m| — g — 1. Hereby it results that the main quantum number 
of these modes, n = rii + n2 — \m\ + 1, accomplishes the selection rules 



|g| < n < \m\ — 1 and k = J\(3\ (^2 — ni + q). 
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6 Concluding remarks 



The study of the Schrodinger equation in Euchdean Taub-NUT space 
is well motivated. The Killing tensors of the Taub-NUT geometry imply 
the existence of extra conserved quantities, quadratic in four-velocities. The 
Taub-NUT case is analogous to the Coulomb problem where an extra de- 
generacy is present and there is the possibility of separating variables in two 
different coordinate systems. This analogy is not perfect since in the Taub- 
NUT geometry the selection rules of the angular quantum numbers depend 
on q and the boundary conditions allow irregular modes. However, this has 
nothing surprising in general relativity where we know already the irregular 
modes of spinless [l^ or spin-half [0 particles in anti-de Sitter space-times. 



In the last time, Iwai and Katayama |T8l 20, pl^] extended the Taub 



NUT metric so that it still admits a Kepler-type symmetry. This class of 
metrics, of course, includes the original Taub-NUT metric. In general the 
Killing tensors involved in the Runge-Lenz vector cannot be expressed as a 
contracted product of Killing- Yano tensors. The only exception is the original 
Taub-NUT space p2[. It will be interested to investigate the Schrodinger 



equation in the more complicated case of the generalized Taub-NUT spaces. 
The extension of the results for this class of space-times will be discussed 
elsewhere [p3|]. 



A 50(3) (g) U{1) harmonics 

The 5*0(3) ® Uil) harmonics defined by Eqs.(|^)-(0) can not be ex- 
pressed in terms of usual spherical harmonics. These are new harmonics 
that must be separately studied by solving the system of eigenvalue equa- 
tions. We start with 

yUO. X) = ^ e^,„(co.^)e^'"^e^^^ (64) 

since then Eqs. (pOf) and (^) are satisfied. It remains to solve only Eq.(|2^) 
and calculate the normalization factor of according to the condition 

= 2 . (65) 



/ d{cos( 



e?^(cos( 
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To this end, we introduce the new variable z = sin^ 9/2 in Eq.(p8D where L 
is given by (|26|). Thus we obtain 



+ (1-2.) A + /(/ + !) 



(66) 



2(1 - 2z)mq 



4^(1 



0. 



This equation has solutions of the form 

e^„(z) ~ z^{l -z)''F{p + k-l,p + k + l + l,2p+l, z) 



(67) 



where the Gauss hypergeometric functions, F, depend on the real parameters 
p and k which satisfy p^ = {m — g)^/4 and fc^ = (m + qY 1^- 

First of all we observe that the solutions (^) are square integrable only 
when F are polynomials selected by a quantization condition since otherwise 
F becomes strongly divergent for z ^ 1. This means that I — p — k must 
be a non-negative integer number. If we replace the functions F of (|67|) 
by Jacobi polynomials [^|, we observe that the solutions of Eq. (|66|) remain 
square integrable if 2p > —1 and 2k > —1. Like in the case of the usual 
spherical harmonics, the good choice is p = (|m| — q)/2 and k = (|m| + q)/2 
where 



1 < Iml < / . 



1^1 

Then by using the normalization condition ( |5B[ ) we find the final result 

{I - \m\)\ {I + \m\)\ 



(68) 



COS! 



V2r+ 

2H 
X (1 — cos 



r(/-g + l)r(/ + g + l) 

m|4-g 



(69) 



2 (l + cos( 



p, 



{\m\-g, \m\+q) 



l-\m\ 



cos( 



For m = \m\ the S0{3) ® f^(l) harmonics are given by ( |64|) and ( |69|) while 
for m < we have to use the obvious formula 



y-Q 



(70) 



When the boundary conditions allow half-integer quantum numbers / and m 
then we say that the functions defined by Eqs. (|6^ and (^) (up to a suitable 
factor) represent SU{2) U{1) harmonics. 
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Thus we have obtained a non-trivial generalization of the spherical har- 
monics of the same kind as the spin- weighted spherical harmonics |2^. In- 



deed, if /, m and q = m' are either integer or half-integer numbers then we 
have 

yi%iO,<P,x) = ^^^^Dl^m'{<P.O,x) (71) 

where -0^^' matrix elements of the irreducible representation of 

weight / of the SU{2) group corresponding to the rotation of Euler angles 
(0, 9, x)- What is new in the case of our harmonics is that these are defined 
for any real number q. 
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